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We present a rigorous formulation of generalized Kohn-Sham density-functional theory. This 

provides a straightforward Kohn-Sham description of many-body systems based not only on particle- 

Qs ' density but also on any other observable. We illustrate the formalism for the case of a particle-density 

based description of a nonrelativistic many-electron system. We obtain a simple diagrammatic 

_,^ , expansion of the exchange-correlation functional in terms of Kohn-Sham single-particle orbitals and 

'"Tt ' energies; develop systematic Kohn-Sham formulation for one electron propagators and many-body 

[•-[- ^ I excitation energies. This work is ideally suited for practical applications and provides a rigorous basis 

for a systematic development of the existing body of first-principles calculations in a controllable 

C^ ' fashion. 
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I. INTRODUCTION 



2 ' Density-functional theorylil allows one to study the ground state properties of the many-body system in terms of 

the expectation value of the particle-density operator. In principle, it offers the possibility of finding the ground state 
energy Eg by minimizing the energy functional that depends on the density only: 



Eg^m:inE[n]. (1) 



^■"^ , Similar theories can be formulaled in terms of the expectation values jOf a spin-density operator or current-density 

, ■ operator, known as spin-densityH and current-density functional theory.B In general, one can imagine a description of 

^ , a many-body system in terms of the expectation value of^aiiy other suitable operator. Such a general description can 

[ — ■ be elegantly presented via the effective action formalism,DO thus leading to a generalized density-functional theory - 

t:]" \ a theory that allows a description of many-body systems in terms of the expectation value of any suitable operator. 

Cn Proving the existence of such a theory is not a trivial matter. In fact, the conventional density-functional theory 

vJH relies heavily on the theorem of Hohenberg and Kohnjd which shows that there exists a unique description of a many- 

J^ ■ body system in terms of the exuectation value of the particle-density operator. Finite-temperature extension of this 

f-^ ' theoreHi|-|Was given by Mermin.El Similar results have been obtained for spin-density and current-density-functional 

■"■^ — , theory.Ha Clearly, it is absolutely important to establish the corresponding existence theorems when constructing a 

j^ ' generalized density-functional theory. In the framework of the effective action formalism, the proof of existence was 

given only in a diagrammatic sense and was tightly lailored to the particular features of the particle-density based 

description of a nonrelativistic many-electron system.Q In this work we present a different resolution of this important 

■ O ■ issue. Our formulation is valid for a general case and does not rely on any perturbative expansions. Not only does 

C , this offer an alternative proof of the previous results for density-, spin-density-, current-density- functional theories, 

but it also provides a rigorous foundation for generalized density-functional theory. 

However important the formal question of existence is, it is of little help for the actual construction of the required 

density or other functionals. For example, to make any practical use of the conventional density-functional theory, an 

k>( I explicit (perhaps approximateLexpression for the energy functional ([I]) is necessary. An important contribution here 

5_] ■ was made by Kohn and ShamE They proposed a certain decomposition of the energy functional, which for a typical 

C^ ' nonrelativistic many-fermion system. 



H = / dx7^t(x) ("-J-V^ +u„„(x) ) ^(x) 

^t(x)^t(y)^(y)^(x) 



dxdy, (2) 



takes the form 



E [n] = Ts [n] + / Uion(x)n (x) dx 

"-W"-(y) 



|x-y| 



dxdy + Exc [11] ■ (3) 



Here Tg [n] is the kinetic energy of an auxiliary system of noninteracting fermions that yields the ground state 
density n (x) , and E^c [n] is the exchange-correlation energy. Once the approximation for E^c [n] has been decided, 
the minimization of the functional (0) leads to the familiar Kohn-Sham single-particle equations.a This approach 
represents the so-called Kohn-Sham density-functional theory. 

A natural question now arises: is there an analog of Kohn-Sham denary-functional theory in the effective action 
formalism? In the original work of Fukuda, Kotani, Suzuki and YokojimaQ the relationship between the two methods 
was not established. Understanding the place of Kohn-Sham theory in the effective action formalism is impogtatit for 
a number of reasons. A vast majority of first-principles calculations are based on the Kohn-Sham method.@'t2l The 
remarkable success of these calculations points to the fact that this seemingly ad hoc decomposition provides a very 
good approximation of the energy functional. As we show in this work, the connection betweeiubhe effective action 
formalism and the Kohn-Sham method can be rigorously established via the inversion method.EiJ The realization of 
this fact immediately leads to a generalized Kohn-Sham theory. 

A distinct feature of the effective action formulation of generalized Kohn-Sham theory is that it provides a sys- 
tematic way of calculating the required functionals. For example, in case of particle-density based description of 
nonrelativistic many-electron systems, it leads to a set of simple diagrammatic rules for constructing the exchange- 
correlation functional entirely in terms of Kohn-Sham derived quantities. We construct the first few orders of the 
exchange-correlation, functional, comment on the local density approximation, and discuss our results as compared to 
other methods .Eiyy 

Applications of the presented formalism is not restricted to the exchange-correlation functional only. We have 
already demonstrated that this method is capable of establishing rigorous Kohn-Sham density-functional formulation 
of one-electron propagators.llJ We briefly discuss the raa.iu_.results of that work and compare it with the existing 
strategies for construction of one-electron propagators .113 tj Lastly, we analyze the excitation energies within the 
effective action formalisro—and comment on the relationship with similar resultstallB obtained via time-dependent 
density-functional theory.L3 

II. EFFECTIVE ACTION FUNCTIONAL 

The effective. action formalism in the context of density-functional theory was discussed by Fukuda, Kotani, Suzuki 
and YokojimaD Here we describe the main features of this method and prove the generalized existence theorems. We 
start by defining the functional W [J] as, 



g-/3iy[J] ^ rp^ (g-0{H+Jil)Q{l)) 



)■ 



Here H denotes the Hamiltonian of the system under consideration. Parameter /3 can be identified with inverse 
temperature. Q is the operator whose expectation value will serve as a main variable of the theory, and J is the 
external source coupled to it. Both H and Q are assumed to be time- independent. Summation over repeated indexes 
is assumed, and the notation J (I) (5(1) embodies all the necessary summations and integrations. For example, to 
formulate the theory in terms of the expectation value of the particle-density operator h (x) , we choose Q — n (x) and 

J(l)(5(l) = I dxJ{x)h{x). 

Thermal expectation value of Q, 

TrfQ(l)e-''(^+''(i')Q(i')) 

Q(i)-^^ 

Trie 
can be written in terms of W [J] as. 



-p{H+J{l')Qil')) 



^(l).^^M (4) 



Denoting the set of allowable external sources as J and the set of all generated expectation values of Q as Q, we can 
establish a map J ^^ Q. We assume that for a given element of J there corresponds only one element of Q, in other 
words, given the external source we can unambiguously establish the expectation value that it generates. Whether 
the converse of this statement is true remains to be proven. The following property of the functional W [J] is of 
fundamental importance. 

Theorem 1: The functional W [J] is strictly concave, i.e. for any a, < a < 1, and J ^ J' 

W [aJ + (1 - a) J'] > aW [J] + {l~a)W [J'] . 

Based on strict concavity of W [J] one can prove the following result. 

Corollary 1: The map J -^ Q is one-to-one. 

The proofs of these statements are given in the appendix. Corollary 1 guarantees that the functional relationship 
between J and Q can be inverted: 

^^ = Q(1)^ J = J[Q]. (5) 

When Q = n (x) , Corollary 1 represents an alternative proof to the theorem of Mermin.Q 

The functional W [J] provides a description of the physical system in terms of the external probe J. We, on the 
other hand, want the description iu— terms of Q. The change of variables from J to Q can be accomplished via a 
functional Legendre transformation.cd This leads to the definition of the effective action functional: 

r[Q] = W[J]~J{l')Q{l'). (6) 

Here, J is assumed to be a functional of Q from Eqn (0) . The functional T [Q] possesses the following important 
property. 

Proposition 1: The effective action functional T [Q] defined on the set Q is strictly convex. 

Proof of this statement is given in the appendix. Differentiating m) with respect to Q, we obtain 

Since our original system is recovered when J — 0, we arrive at the important variational principle: the functional 
r [Q] reaches a minimum at the exact expectation value of Q, 

In the zero temperature limit, (3 — > oo, Qg represents the exact ground state expectation value of Q and T [Qg] equals 
the exact ground state energy. Obviously, at finite temperatures 

T[Qg] = -^\nTr(erP"), 



ig 



(1) 



Tr(g(l)e-'^^ 



Tr e-'9« 



The effective action formalism furnishes a rigorous formulation of generalized density-functional theory. The exis- 
tence of r [Q] is guaranteed by Corollary 1. When Q is a particle-density operator, we obtain conventional density- 
functional theory; when Q stands for a spin-density operator-, we have spin-density- functional theory ;EI when Q is a 
current operator, we obtain current-density-functional theory,l3, etc. 

To implement this formally exact method, an approximation of the effective action functional T [Q] is required. 
For a particle-density based description of a nonr|elativistic many-electron system the efEective action functional can 
be approximated via the auxiliary field method.u However, this method does not leadu to Kohn-Sham version of 
density-functional theory. The relationship between_the effective action formalism and Kohn-Sham density-functional 
theory can be established via the inversion methodJlil 



III. GENERALIZED KOHN-SHAM THEORY VIA THE INVERSION METHOD 

Consider the following general Haniiltonian 

H = Hq + XHint 

which depends on the coupling constant A as a parameter. The same is true for the effective action functional 

r = r[Q,A]. 

Clearly Q, A are to be considered as two independent variables. Note, however, that this does not prevent the exact 
expectation value Qg from depending on A: this dependence is fixed by the variational principle 

The functional T [Q, A] is defined as, 

r[Q,A] = M^[J,A]-J(l')Q(l'), (9) 

where J is functional of Q and A. This functional dependence is provided by the equation 

5 J {I) ^^^- 
The inversion methocH'tll proceeds by expanding all the quantities in Eqn (S) in terms of A; 

J [g. A] = Jo [Q] + A Ji [0] + AV2 [Q] + . . . , 

W [J, A] = Wo [J] + XWi [J] + \^W2 [J] + . . . , (10) 

r [Q, A] = To [Q] + AFi [Q] + A^Fa [Q] + . . . . 
Comparison of the two sides in Eqn (g) for different orders of A, 

5^AT,[Q] -^AW. [^A^Jfe[Q]" 

-5]AV.(1)Q(1), (11) 

leads to the expression for F; [Q] , 

Ti [Q] = Wi [Jo] + J2 ^-%T^J>^ (1) - J^ (1) Q (1) 



I ki + ...+k, 

y — y 



I _, fci + ...+fc„<i rmil/ r T 1 

^l-{kx + ... + k^) [JoJ 

X 



TO! •' — ' (5Jn (1) . . . (5 Jo ("^) 

xjfc^ (I)--- J/c„ (rn) 

Functionals {Wi [Jo]} and its derivatives are assumed to be known; they can usually be obtained via standard many- 
body perturbation techniques (specific examples will be given in the next section). Since Q and A are considered to 
be independent, it follows from Eqn (R) that functionals {J^ [Q]} can be obtained using. 

Consider the zeroth order term, 

ro[Q] = Wo[Jo]-Jo(i)Q(i). (13) 



Using Eqn (|T|) 



SWo[Jo]SJ,{V) nn'^^Lill 



fswoiJo] Q^r))'-^^o. 



\6Jo{V) ^' 7 6Qil) 

Strict convexity of Fg [Q] (see Proposition 1) prohibits {SJq (!') /6Q (1)) from having zero eigenvalues. Thus we obtain 
that Jo obeys the equation: 

«(')^^- <») 

Hence Jq is determined as a potential which generates the expectation v&lue Q in the noninteracting (A — 0) system. 
Notice that the same exact notion appears in Kohn-Sham formalism.EI We refer to this noninteracting system as 
Kohn-Sham (KS) system and Jq as Kohn-Sham potential. Eqn (111) allows one to simplify the expression for F; [Q], 
which now becomes, 

r^ [Q] = Wi [Jo] - Si.oJo (1) g(l) + E ^-^TT^Jk (1) 



rn=2"^' fc„.'X.>l '5^o(l)...<5Jo(m) 

xjfcj (I)--- Jfc„ (m-). 






The important message here is that the expression for F; [Q] involves only lower order functionals (Z — 1, Z — 2, . . . , 0). 
Thus starting with Jo we can determine Fi [Q] as 



From Fi [Q] we can find Ji as, 



Ji(l) 



Fi [Q] = Wi [Jo] . (15) 

_aMO] ^ (5Fi [Q] (5Jo {!') 
6Q (1) 5Jo (1') 6Q (1) ' 



^Jo (2) ' 
where the inverse propagator is defined as, 

5 Jo (2) _ ( S'^Wo [Jo 



^i(l)-^(l,2)?^, (16) 



P(l,2) = - 



^Q(l) V'^Jo(l)<5Jo(2) 

Once Ji is known, we can find F2 [Q] : 

dJo (1) 



or using Eqn (|lq ) 



2Wo(l)Wo(2)^^^^)'^^(2)' 



-l->4^-''.^)tf' 



From r2 [Q] the expression for J2 follows as, 

SW2 [Jo] 



J2(l) -^(1,2) 



SJo (2) 



This, in turn, leads to Fs [Q] and so on. In this hierarchical fashion one can consistently determine F [Q, A] to any 
required order; 

To ^ ^0 ^ Ti ^ ,/i ^ F2 ^ J2 ^ • • ■ ■ 

The important point here is that all higher orders are completely determined by the Kohn-Sham potential Jq. Let us 
now apply the variational principle to our expansion of F [Q, A] : 

sr[Q,x] _Q 



Since 



SQil 
STo [Q] 



-^o(l) 



SQil) 
we have the following important result 

J /IN _ ST,nt [Q] _ , M SV„U [Q] . „^ 



where 



^^nt [Q] - Y. ^'r.: [Q] ■ 



Since at any order the effective action functional is completely determined by Kohn-Sham potential Jo, we arrive at 
generalized Kohn-Sham self-consistent method: 

1. Choose the approximation for Tint [Q] (one obvious choice is to truncate the expansion at some order). 

2. Start with some reasonable guess for the Kohn-Sham potential Jo- 

3. Calculate Tint [Q] ■ 

4. Determine new Kohn-Sham potential Jo via Eqn (M) . 

5. Repeat from step 3 until self-consistency is achieved. 

The formalism described above can be applied to any general case and provides a rigorous basis for generalized 
Kohn-Sham theory. Practical implementation of the self-consistent procedure obviously depends on the particular 
Hamiltonian under consideration and the choice of the operator Q. In the next section we apply this method to the 
case of a particle-density based description of nonrelativistic many-electron system. 

IV. KOHN-SHAM DENSITY-FUNCTIONAL THEORY. 

A. Derivation of Kohn-Sham decomposition 

Let us consider a typical nonrelativistic many-electron system described by the Hamiltonian (0) and develop the 
description in terms of the particle-density operator: 



Q(l) = ^t(x)^3(x) = n(x). 



For convenience, the spin degrees of freedom are suppressed here (those can be easily recovered if necessary). The role 
of the Gcuipling constant A is played by e^. We now evaluate the effective action functional F [n] using the inversion 
methodO described in the previous section. Coupling constant expansion of the functional W [Jo,e^] , 

W [Jo, e'] = Wo [Jo] + e^Wi [Jo] + e^W2 [Jo] + ■ ■ ■ , 
can be conveniently generated using path integral representationE3 



Here 



S[ilj'<,i}] = / dxi^^x) 



dr 2r 



'^hon \P^) 



Tp{x) 



dxdx' ^p^ {x)tlj' {x')u{x ~ x')^p{x')'tp{x), 



(20) 



(21) 



u{x — x') 



5{t- t') 



and V'^jV' denote Grassmann fields.E3 
The zeroth order term is given by, 



dx = I dr I dx. 



Wo[Jo]=-^J2^n{l + e-^'^), 



where e^'s denote single-particle energies of Kohn-Sham noninteracting system: 

V2 



--Z h Vioni^) + Jo (x) ) (fii (x) = Eiipt (x) 



*W = ^f^^\f^(.^)\'^ ^ 



n,: = fe^-' 



+ 1)" 



First order term is given by, 



VFi [Jo] 



Here solid lines denote Matsubara Green's function of Kohn-Sham noninteracting system 



5o(xr, x'r') 



T > r 



,E¥'»(x)¥'*(x')e-^'("""')(«.-l) 



(22) 



and dotted line stands for the Coulomb interaction u{x — x'). The explicit expression for Wi [Jo] is given by. 



Wi[Jq] ^-^11 Go{x,x)u{x - x')Qo{x',x')dxdx' 



"na I f Qo{x,x')u{x ~ x')Qo{x\x)dxdx' . 



(23) 



The second order term is 



W2 [Jo] = 




(24) 



(25) 



(26) 



Similarly, one can generate higher order terms. Let us now calculate the first few orders of the effective action 
functional. Using Eqn ([13) , the zeroth order correction is given by, 



-Pei^ 



In the zero temperature limit, /? ^ cxd, this transforms to 

N 



Jq (x) n (x) dx. 



^0 [n] = X! ^* ~ / Jo (x) n (x) dx, 



"0 [n] = To [n] + ', 



i(x)n(x)dx. 



Hence at zeroth order, the effective action functional is given by the sum of kinetic energy To [n] and ion-potential 
energy of the Kohn-Sham noninteracting system. From Eqn dTS) , the first order correction is given by, 

Ti [n] = Wi [Jo] . 

The expression for Wi [Jq] has been presented earlier (see Eqn (|2^)). In the zero temperature limit 



El \n] = - 
iL J 2 



n(x)n(x') 



dyidx! 



x-x'l 

n (x,x') n (x',x) 



dxdx'. 



where n(x,x') = ^^i (x) (^* (x') . Here, the first term represents the classical Hartree energy and the second term 
represents the exchange energy. Both are evaluated with respect to Kohn-Sham noninteracting system. 

Postponing the evaluation of the second order term until the next section, let us summarize the results we have 
obtained. It is clear that in the zero temperature limit the expansion for the effective action functional. 



F [n] = To [n] + / Uion(x)n (x) dx 



n (x) n (x') 



dxdx' 



n(xxOn(x',x) ^^^^, 



i=2 



e'T, [n] 



(27) 



coincides with the deconippsition proposed by Kohn and ShamE Therefore the inversion method of evaluating the 
effective action functionallill naturally leads to the Kohn-Sham density- functional theory. Application of the variational 
principle (see Eqn (0) ) to the expansion (p^) yields the well-known Kohn-Sham self-consistent procedure and the 
corresponding single-particle equations. Comparison of Eqns ^ and (|2j) immediately provides an expression for the 
exchange-correlation functional. These topics are discussed in detail in the following sections. 

To conclude this section, weuvould like to note that alternatively, the effective action functional can be also evaluated 
via the auxiliary field methodu However, in this method thje Kohn-Sham decomposition has to be artificially imposed 
to allow the study of the exchange-correlation functional or Kohn-Sham density-functional theory in general. 



B. Construction of the exchange-correlation functional 

The success of first-principles calculations based on Kohn-Sham density- functional theory depends on the accuracy 
of the approximations to the exchange-correlatimi functional. The analysis of Kohn-Sham theory, or E^c in particular, 
via standard many-body perturbation theoryQ'E3 was always a challenging task, for there was no explicit connection 
between the two methods. The advantage of the effective action formalism is that it is a rigorous many-body ap- 
proach specifically designed for a density-based description of many-body systems. This formalism provides a natural 
definition of the exchange-correlation functional as. 



This expression involves only Kohn-Sham based quantities and is especially suitable for practical appiieations. Simple 
diagrammatic rules for evaluating higher order terms in the expansion of T [n] are readily available.uO This, in turn 
leads to the following set of rules for the calculation of the exchange-correlation functional: 

1. Draw all connected diagrams made of Kohn-Sham propagators Qq {x, x') and Coulomb interaction lines u{x — x') 
with the corresponding weight factors Bj 

2. Eliminate all the graphs that can be separated by cutting a single Coulomb interaction line. 

3. For each two-particle reducible (2PR) jgraph (i.e. any graph that can be separated by cutting two propagator 
lines) perform the following procedures 

(a) Separate the graph by cutting 2PR propagators. 

(b) For each of the two resulting graphs join two external propagators. 

(c) Connect the two graphs via the inverse density propagator P (x,x') . 

(d) Repeat the procedure until no new graph is produced. 

(e) Sum up all the resulting graphs including the original graph. 

The inverse density propagator 2? (x, x') is given by. 



2?(x,x') 
or in terms of Kohn-Sham orbitals 

2?(x,x') 

Application of these rules for gq 



/ go(xT,xV)go(xV,XT)dT' 
-'0 






LP, (X) if* (X') V?j (X') LP* (x) 



(28) 



Kohn-Sham exchange functional :E5EjE^ 



-of the first-order correction to Exc [r^ is obvious and leads to the well-known 



E. 



XC.l 



-m 



Let us now consider the second order correction to Exc [n] ■ Rule # 1 leads to the already given expression for W2 [Jo] 
(see Eqn (Eq)). Application of Rule # 2 leads to elimination of the last two graphs, thus giving 



The third graph in the above expression is 2PR. According to Rule # 3 it transforms to: 




Here double solid line denotes the inverse density propagator ^(x, x'). Therefore the final expression for the second 
order correction to E^c is given by, 



E.rr_2 — -tt: 




It is instructive to apply the above procedure to the case of a homogeneous electron gas in the zero temperature limit. 
It can be demonstrated that in this case Rule #3 leads to the complete elimination of 2PR graphs. Indeed, at this 
limit one can show that 



*> ' *CrO 



= 0. 



Here, black circles denote parts of the diagram that are connected to each other via two propagators. For example, 
in the expression for the second order correction to E^c the last two graphs completely cancel each other 



and as expected 



e 




jTihom ^ 




Local Density Approximation (LDA) represents a popular choice for Exc [n] in first*rinciples calculations. This 
approximation and subsequent corrections can obtained via the derivative expansioncj of the exchange-correlation 
functional 



Ex 



ii;(";(n(x))+42)(n(x))(Vn(x))^ 



dx 



where Exc is a function of n (x) , not a functional. 
LDA corresponds to the first term in this expansion; 



Si-f^M^ EtHn{^))d^. 



Function ExJ {n (x)) can be found by evaluating the above expansion at constant density n (x) — no 



Ei°J im) = yExc [no] 



no = const. 



Here Exc [no] represents the exchange-correlation energy of the homogeneous electron gas with density no- 

Regarding the results obtained in this section, we would like to emphasize the following points. Not only does the 
effective action formalism leads to a straightforward set of rules to calculate the exchange-correlation functional up 
to any arbitrary order, but it can also be used to generate similar quantities for descriptions based on the observables 
other than particle density (for example, current-density functional theory). One should also remember that the 
expansion of Exc [n] represents only part of the general picture provided by the effective action formalism. As we 
show later, the same formalism allows us to develop a rigorous and systematic Kohn-Sham theory for one-electron 
propagators and many-body excitation energies. 



10 



C. Kohn-Sham self-consistent procedure 

As we have demonstrated earlier, application of the variational principle leads to the Kohn-Sham self-consistent 
procedure. In the case of the traditional density-functional theory, the typical self-consistent procedure takes the form 

1. Start with some reasonable guess for the Kohn-Sham potential Jq (x). 

2. Solve Kohn-Sham single-particle equations. 

3. Determine new Kohn-Sham potential Jq (x) using 

Jo (x) = / ,^ L d^' + ^xc (x) , 

where the exchange-correlation potential Vxc (x) is defined as, 

5E,;c [n] 



Vxc (x) = 



4. Repeat from step 2 until self-consistency is achieved. 



5n (x) 



In LDA the exchange-correlation functional is an explicit functional of electron density n (x) , and the exchange- 
correlation potential Vxc (x) can be obtained by straightforward differentiation. In the case of the diagrammatic 
expansion, this simple property no longer holds and Exc [n] appears as an implicit functional of n (x) . The exchange- 
correlation potential Vxc (x) can still be found using the equation: 



„„(X) = -/P(X.X',^^'. 



The functional derivative in the above expression can be easily evaluated based on the following relationships: 

5Qq {xi,X2) f'^ 



oJa (x) Jo 



" ^ ' > ^ 22? (x,xi) 00(2^1, 0:2) P(X2,X'). 




SGo (a;i,a;2 
For example, the first order correction to Vxc (x) is given by, 

VxcA (x) = - 

Obviously, even in case of the diagrammatic expansion of Exc one could still use the above mentioned self-consistent 
procedure. To reduce the computational effort, however, slight modification of that procedure might be advantageous. 
Namely, we suggest to shift the emphasis from density n to Kohn-Sham potential Jq- Indeed, Corollary 1 guarantees 
that there is a one-to-one correspondence between n and Jq. Thus, we can consider the effective action functional 
that depends on Jq rather than n: 

r [Jo] = r[n [Jo]] . 

The variational principle then takes the form 

^r[Jo] 



SJo (x) 



= 0. 



In other words, one has to find a Kohn-Sham potential that minimizes the effective action functional r[Jo]. To 
accomplish this task, one could use the so-called steepest descent minimization method.EZI In this case the self- 
consistent procedure takes the form: 

11 



1. Start with some reasonable guess for the Kohn-Sham potential Jo (x). 

2. Calculate the direction of steepest descent s (x) as, 

^'~ Wo(x) 
or 

j\ J,,' ^r.i„f [ Jq] 

SJo (x) 



s (x) = - / P-I (x, x') Jo (x') dx' - 



3. Determine new Kohn-Sham potential Jq"^'" (x) from the old Kohn-Sham potential Jq''' (x) by stepping along 
the direction of steepest descent: 

Jo"^-(x) = J°''^(x)+as(x) 

where a is the length of the step. 

4. Repeat from step 2 until self-consistency is achieved. 

The advantage of this self-consistent procedure is that it is much easier to calculate 

^r[Jo] 

<5 Jo (x) 
rather than 

(Sr \n] 



Sn (x) 
For example, when Tint [Jo] is approximated by its first order correction 



1/ A A X im 





^'"' ["^o] - 2/3 \^ \^ 2/3 ^ly 
the steepest descent direction is given by 

.(x) --Jv-' (x,x') (^Jo (x') - J ^^j^dx") dx' 

The advantage of this expression as compared to the first order correction to Vxc.i is that here, we avoid the calculation 
of the inverse density propagator V (x, x'). Note that V^^ (x, x') can easily be written in terms of Kohn-Sham single- 
particle orbitals and energies (see Eqn (Eq)). 

V. TIME-DEPENDENT PROBE 

To study excitation energies and one-electron propagators, it is necessary to consider an imaginary-time-dependent 
probe. The definition of the functional W [J] is changed correspondingly 



Note that the parameter /3 has been absorbed into W [J]. In order to proceed with the inversion method, we need to 
assure that the map J (x) — > n (x) , 
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is invertible. Since at the end the time-dependent probe is set to zero, we can assume our source to be infinitesimally 
small. Therefore it suffices to prove the invertibilty in the small neighborhood of time-independent external source: 

J (x) + 5 J [x) -^ n (x) + 5n {x) . 

Since, 

sr2T 



we need to show that the operator 



JJ{x)5J{x')J j^^^ 

has no zero eigenvalues. This property follows from the following theorem. 

Theorem 2: The operator W^"^' ix,x') is strictly negative definite. (The proof is given in the appendix.) 
This theorem guarantees that as long as two infinitesimally small probes differ by more than a pure time-dependent 
function, they would produce two different densities. Once this one-to-one correspondence has been established, 
further analysis proceeds similar to the time-independent case. The effective action functional is defined as, 



T[n]^W[J]~ I J (x) n (x) 



dx. 



where J is assumed to be a functional of n by Eqn ( |3(|) . Using the inversion methodjl^ the effective action functional 
can be found as a power series in terms of the coupling constant: 



r [n, e^] = To [n] + e^Fi [n] + e^T2 [n] 



To [n] = Wo [Jo] - / Jo {x) n {x) dx 



Zeroth order term is given by. 



The functional Wq [Jq] describes a Kohn-Sham system of noninteracting electrons in the presence of an imaginary- 
time-dependent external potential Jq {x) . Kohn-Sham potential Jq {x) is chosen such that the time-dependent density 
n {x) is reproduced. Diagrammatic structure of Ti [n] is the same as its time-independent counterparts with the only 
difference being that Kohn-Sham propagator is now defined in the presence of a time-dependent Kohn-Sham potential 
Jo {x) , 

/ 5 V^ \ 

Qq^ {x,x') ^ -[- h Wion(x) -f Jo {x) 5{x - x'), 

\ OT 2m / 

and the inverse density propagator becomes 

v{x,x')^-[goix,x')go{x\x)r\ 



VI. ONE-ELECTRON PROPAGATORS 

Using the effective action formalism it is possible to develop a systematic Kohn-Sham density-functional approach 
to one-electron propagators, ft is common practice to use converged Kohn-Sham single-particle orbitals and energies 
in quasiparticle calculations. However, very often these methods are not very systematic in their use of Kohn-Sham 
based quantities.EZI The formalism presented below provides a rigorous theoretical foundation for the calculation of 
quasiparticle properties based on Kohn-Sham noninteracting system. 

Consider W [J] in the presence of the auxiliary nonlocal source ^ (x, x') 

f3 



where 

The nonlocal source, ^(x,x'), allows us to write the one-electron propagator or the (finite temperature) Green's 
function, Q{x^ x') = — (Tr'4' {x) V'^ {x')), as a functional derivative 

Using the well-known property of the Legendre transformationJHj 

5W [J]\ _( 5T [n] \ ^g^^ 



^5i{x',x)Jj \5^{x',x)^ 
the one-electron propagator can be expressed in terms of the effective action functional as, 

«(.,.')^(j^) +(ff^) . (32) 

Using the property (Bl]) we obtain that 



xci I \ I \ xcf > \ I =Goix,x'). (33) 

'5Ua;,a;)/„ \dC{x',x)J j^^ 



Let us consider the second term in Eqn (|32|). One can show that 



' - ' ' -go{y,x)go{x,y)dydy 



\SC{x',x)J ^^ J J Sgoiy,y 

STmt [n 
6n{y) 



Qo{y,x')Qti{x,y)dy. (34) 



Our original system is recovered by setting i^ to zero. Using Eqns (|32|), (|33|), (|34| ) the exact one-electron propagator 
(in operator notation) is given by, 

Q = Qo + Qo-^o- Go, (35) 

where the Kohn-Sham self-energy Eq is given by, 

So(xi,X2) = ^^ '"* - Jo{xi)5{xi - X2) (36) 

oyo (a;2,a;i) 

and Jo(a;i) is the Kohn-Sham potential, 

ST^nt [n] 



Jo{x) = 



Sn{x) 



The functional derivative in Eqn (36) can be easily evaluated since the functional Tint can be expressed entirely in 
terms of Kohn-Sham Green's functions Qq. The above expression for the one-electron propagator can also be written 
as, 

G = Go + Go ■'^0 ■ G, 
where 
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Eo = So ■ (1 + 00 • So)-' = i^o' +5o) ' . 
Therefore the exact self-energy S is given by, 

S (x, a;') = Jo{x)S{x — a;') + ^a{x, x'). 

The above formulation provides a systematic way to study one-electron propagators and self-energy in Kohn-Sham 
density- functional theory. The important feature of this formulation is that both self-energy and Kohn-Sham potential 
arc determined from one quantity Tint [n] ■ A single approximation to the functional Tint [n] simultaneously generates 
both the self-energy and Kohn-Sham potential. This is to be contrasted with the common strategy of performing 
separate approximations for Kohn-Sham potential and the self-energy.E3 

VII. EXCITATION ENERGIES 

In addition to one-electron propagators, pthe effective action formalism allows a systematic study of the many-body 
excited states in density-functional theory.Q Consider the Fourier transform of M^(^) (x, x') : 

Ty(2)(x,x',zi/,) ^ / W'-''^{x,x')e'''^V-^)d(T-T'), 
Jo 
Vs = 27rs//3, 

It can be analytically continued into the complex plane u: 

VF(2) (x,x',cj) - W^(2) (x,x',z:/,) |„,^^+.^. 

The proposition below guarantees that this analytic continuation has an inverse in the upper complex plane lj including 
the real axis. 

Proposition 2: The operator VF^^^(x, x',a;) has no zero eigenvalues when uj is located in the upper half of the 
complex plane including the real axis, (see appendix for the proof.) 

Let us define the excitation kernel as, 

Ti^H^,x')-f ^^^[^] 



^(5n(x)^n(x')/„(^) 
It easy to show that 

fw(^\x,x')T<^^\x\y)dx' = -Six-y), 

or in terms of Fourier transforms 

VK(2) (x, x', ij.,)r(2) (x', y, -iiys)d^' = -5 (x - y) . (37) 

Proposition 3: r^^^(x, y, —ii^s) has a unique analytic continuation r'^^^(x, y, — w) such that 

1. it does not have any zeros in the upper half of the complex plane lo including the real axis. 

2. [r(2)(x,y,-c.)]-'^^, (c^^). 

The proof is given in the appendix. Based on the above statement we can extend the relationship (B^ to the whole 
complex plane u> : 

J W^(2) (x^ x'^ ^)r(2) (x', y, -a;)d^' = -5 (x - y) . (38) 

Consider VF'^^(x, x',a;) in the zero temperature limit (/3 -^ cxd). It is well-known that it has poles just below the real 
axis at the exact excitation energies, 
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UJ = El - Eq - IT]. 

Then it follows from Eqn ( pSj ) that r(^^(x, y, — w) has zero, eigenvalues at the exact excitation energiesJfl In other 
words, when lu = Ei — Eq — irj , there exists ^ (x) such that,Q 



yr(2)(x,y,-c.)^(x)dx = 0. 
Using the coupling constant expansion for T(n) the excitation kernel becomes 

Substituting this in Eqn ( p9|) and using the fact that, 



(39) 



r[,')(x,x',-c.) 



.(2) 



w^'>{^y,uj), 



we obtain 



r{2) 



e(x)= / /W^f(x,x',c.)rl^i(x',y',-c^)ay')rfx'dy' 



(40) 



Here Wq (x, x'jCj) represents the negative of the density-correlation function for Kohn-Sham noninteracting system; 



W^o(')(x,x',^) = ^ 



nl (x) ni (xQ _ ni (x) n* (xQ 
u) — ujf'^ + ir] uj + ujf'^ + irj ' 



where 



and 



(Z'=^|n(x)|0^^)=n,(x), 



, ,ks jT'ks TT^ks 

LOi ^ El - Ef) . 

Here all the quantities refer to the Kohn-Sham noninteracting system. Aside, from the notational differences, Eqn 
( [40|) coincides with similar expressions derived using time-dependent DFT£a'E3 Searching for a solution of the form 

^ (x) = ^ ami (x) + Y^ hiril (x) , 

we obtain the following matrix equationjl3 

n/f 1 r /I 1 

= (cj + i-q) 



L M 
M* L* 



A 
B 



-I 
1 



A 
B 



where 



M,j^ 



I J n* (x) ri^l (x, y, -u;)n* (y) dxdy, 



^i (x) r,Vt(x, y, -i^)nj (y) dxdy 



(2) 
{Ei ~ Eq) Sij, 



and 



The values of lu for which the above matrix equation has nontrivial solutions, is determined by the condition 



det 



L M 
M* L* 



{uj + irj) 



-1 
1 



= 0. 



Therefore, the effective action formalism presents an alternative way (as compared to time-dependent density- 
functional theory) for calculating the excitation energies. However, in addition to it, the effective action formal- 
ism also provides a means of calculating|-t|he exchange-correlation kernel. This feature is missing in conventional 
time-dependent density-functional theory.E3 
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VIII. CONCLUDING REMARKS. 

Based on the effective action formalismQ we developed a rigorous formulation of generalized Kohn-Sham theory. 
This formulation is specifically geared towards practical calculations of the ground and excited properties of real 
systems. Indeed, in the case of a particle-density based description of nonrelativistic many-electron system we arrived 
at a systematic way to study the exchange-correlation functional, one-electron propagators and many-body excitation 
energies entirely in terms Kohn-Sham single-particle orbitals and energies. The presented formalism is very general 
and can be applied to a various many-body systems for constructing Kohn-Sham like description in terms of the 
expectation value of any general operator (e.g. spin-density, current-density- functional theory). 
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APPENDIX A: PROPERTIES OF FUNCTIONALS W [J] AND F [Q] . 
1. Time-independent probe 

Theorem 1: The functional W [J] ~ —\ InTrf e^'^v'^^'^'^^'^^^'^ J is strictly concave, i.e. for any a, < a < 1, and 

W [a J +(!-«) J'] > aW [J] + il-a)W [f] 

ProofS 

Consider two Hermitian operators A and 13. Let {^i} be a complete set of eigenstates of the operator aA+{l — a) B. 
Then 



Tr (e"^+(i-")^) =^(*. 



.^aA+{l-a)B 



*, 



y^ ga(*,|i|*,) + (l-a)(*.|B|*.) 



Using Holder's inequalityEj (a^, 6i > 0) , 



^a?r"< E-0 E^^ 



we obtain that 



Tr (e"^+(i-")^) < |^e(*-l^l*') ) (^e(*'l^l*' 



(l-a) 



From Holder's inequality it follows that the equality in the above expression holds only if for any i 



^, 



A 



*,:> = ^ 



B 



*0 + x, 



where x is constant independent of i. 
Since e^ is a convex function, 



(Al) 



i i 


e^ 


* 


y^e(*'l^l*0<^(vl/. 


e^ 


^ 



(A2) 
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Equal sign in the above equations holds only if {$i} are the eigenstates A and B . Collecting all the results together 
we obtain the following inequality: 

Tr (e"^+(i-")^) < (^Tre^) " (Xre^) ^'""^ . (A3) 

Equality holds here if and only if operators A and B differ by a constant: 

A^ B + const. (A4) 

Indeed: 



) Suppose Eqn A4 is true, then equal sign in (A3) is obvious 



) If there is an equal sign in (A3) then operators A and B must have a common set of eigenstates {vti} and 



condition (Al) must be true. Then considering the representation of operators A and 13 in the basis of eigenfunctions 

A = B + const. 



{^i} we obtain 



Therefore 



when A =f^ B + const. 
Setting 



and 



we can easily obtain 



Tr (e"^+(i-")^) < (Tre^)"(Tre^) 

A = -(3(^H + J{l)Qil)) 
B = -I3(h + J'{1)Q{1)), 



(l-a) 



W [aJ + (1 - a) J'] > aW [J] + il-a)W [J'] , 



when J ^ J'. QED. 

Corollary 1: The map J ^ Q, is one-to-one. 
Proof: Consider the functional 

A[j] = w^[j]-j(i')g(i'), 

where J and Q are considered to be independent. Since W [J] is strictly concave (Theorem 1) , it follows that A [J] 
is strictly concave. Therefore if A [J] has an extremum it is unique. Hence if the equation 

U{\) U{\) ^^' 

has a solution, it is unique. QED. 

Proposition 1: The effective action functional F [Q] defined on the set Q is strictly convex. 
Proof: Consider the family of the functional {A [J, Q] : J E J , Q E Q} , 

A[J,Q] = W[J]-J{l')Qir). 

Here Q is considered to be independent of J. Obviously A [J, Q] is linear in Q. The effective action functional can be 
defined as 

r [Q] = sup {A [J, Q],JeJ}. (A5) 

The above expression represe nts the most general way to define the effective action functional. In our case two 
definitions, Eqn (a) and Eqn (A5) , arc equivalent. Convexity of F [Q] follows from the fact that it is a supremum of 



the family of the linear functionals (in Q). Because by construction for any element J E J , there corresponds only 
one Q G Q , the functional F [Q] is also strictly convex. QED. 
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2. Time-dependent probe 



Theorem 2: The operator W^'^' {x, x') is strictly negative definite 

Proof: We have to prove that there exists no function / (x, t) 7^ r] (r) such that 



/* {x) W^^> {x, x') f {x') dxdx' > 0. 



(A6) 



The operator VK^^^ {x, x') actually represents the negative of the density correlation function in the presence of the 
time-independent source J(x): 



VF(2)(x,x') 



/ DiP^Dj^n {x) n [x') g-^-^^^-^l 



where 



[V'^v] -5'[?/'^v] + f Jix)n{x)dx, 



and n (x) denotes density fluctuation operator, 

n {x) = V' {x) ip {x) ~ n (x) . 
Defining a Fourier transform as, 

/■/9 , ,, 

M/^(2)(x,x',it^,) = / W^^'^{x,x')e"''^^-^)d{T-T'), 



Vs = 27rs//3, 



left hand side of (A6) transforms into 

Y. f f(f (^' ''s))* W^^'' (x, x', ti^s) f (x', i/,) dxdx'. 



(A7) 



Strict concavity of W [J] for time-independent case guarantees that 1^3 = term is strictly negative definite. For 
convenience, we ignore this term from the sum. Lehman representation for W^'^' (x, x', Ws) is given by, 



W^^> {x,x\ii^. 






rnl 

X (m \n (x)| I) {I \n (x')| m) , 



E,n - Ei+iVs 



(A8) 



where 



H + / J (x) n (x) \m) = Er,i\m) 



H+ J{^)n{^)] \l)^Ei\l) 



Therefore an arbitrary [vg ^ 0) term in the sum (A7) can be written as. 



^w[j] v^ e ^ " - e *" ' 



E 



Em- Ei+ iVs 



n(x)/(x, i^s)rfx 



.W[J] 



E 

E^<Ei 



n(x)/(x, j/5)dx 



Sm - El+ iVs El - Em + il^s 
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After a simple algebra we obtain 



,W[J] 



E 2 



n (x) / (x,j/s) dx I 

{El — Em) ■ 



g-/3E, _ ^-m. 



{El - E^n) -f 

Obviously, the above expression cannot be positive. Moreover, it cannot be zero since this would imply that for all 
m ^ I 



j n(x)/(x,i/s)dx 



= 



and J n (x) / (xji^s) dx commutes with the Hamiltonian. This is impossible unless / (xji/^) is constant independent 
of X or / (x,t) is a function of t only. This case, however, was excluded from the very beginning. Therefore, the 
operator VF^^) {x^x') is strictly negative definite. QED. 

Proposition 2: The operator VI^"-'(x,x',lj) ~ W^'^' (x, x', w.,) \w^^uj+i7-i has no zero eigenvalues when oj is located 
in the upper half of the complex plane including the real axis. 

Proof: Suppose / (x) is an eigenvector corresponding to a zero eigenvalue then 

(A9) 



/* (x) Ty(2) (x, x', Lu)f (x') dxdx' = 0. 
Using Lehman representation we obtain 

0= E 



E„,<Ei 



n (x) / (x) dx 



g-/3B„ _ g-/3Ei 



-PE„ 



-PEi 



z-{Ei- E^) z + {Ei- E„ 
where z = co + irj. Consider the imaginary part of the above expression. 



= -ARe (z) Im (z) y^ (m n (x) / (x) dx 



E„^<E, 



{El ~ Em) {e-^^- - e-P^') 



z2 - {El - Emf 



Re{z)Im{z) = 0. 



By Theorem 2, the real part of z cannot be zero, since in this case VF'^2^(x, x',aj = ivs) is strictly negative definite 
and condition ( |A9| ) cannot be satisfied. Therefore we necessarily obtain that Iin {z) = and 



Lo + iri = real 



uj ~ real —ir\ 



Therefore ^^''^'(x, x',a;) may have zero eigenvalues only when uj is located in the lower half of the complex plane 
excluding the real axis. QED. 

Proposition 3: r(2^(x, y, — Ws) has a unique analytic continuation r^2)(x, y,a;) such that 

1. it docs not have any zeros in the upper half of the complex plane oj including the real axis 

2. [r(2)(x,y,c.)]~'^^, (^^^). 

Proof: 

Existence follows from the fact that r(2^(x, y,a;) can be defined as — [VF'2^(x, y,a;)] . Suppose there exist two 
analytic continuations of r'^2^(x, y, —ivs) with the properties 1 and 2. Then VF'^2^(x, x', w^) will have two different 
analytic continuations which is impossible. QED. 
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